Group Theory Approach to Band Structure: Scarf and Lame Hamiltonians 



Hui Li and Dimitri Kusnezov * 
Center for Theoretical Physics, Sloane Physics Laboratory, Yale University, New Haven, CT 06520-8120 

(November 30, 1998) 



OS 
OS 
OS 



3 
^1" 



> 

(N 
O 
(N 

o 

Os 
Os 



c 

o 
o 



13 



The group theoretical treatment of bound and scattering 
state problems is extended to include band structure. We 
show that one can realize Hamiltonians with periodic poten- 
tials as dynamical symmetries, where representation theory 
provides analytic solutions, or which can be treated with more 
general spectrum generating algebraic methods. We find dy- 
namical symmetries for which we derive the transfer matri- 
ces and dispersion relations. Both compact and non-compact 
groups are found to play a role. 



The application of group theoretical techniques to 
physical problems has a long and fruitful history. From 
the early works on the low- lying spectrum of hadrons jl| , 
to the more recent studies of nuclear and molecular 
H many-body problems, the theory of Lie groups and 
algebras has become important not only in explaining 
the behavior of various physical systems, but also in con- 
structing new physical theories. Problems of interest can 
usually be approached in this manner when a spectrum 
generating algebra (SGA) can be identified. A SGA ex- 
ists when the Hamiltonian H can be expressed in terms of 
generators of an algebra. As a consequence, the solution 
of the Schrodinger equation then becomes an algebraic 
problem, which can be attacked using the tools of group 
theory. In addition to providing a general approach to 
solving the problem of interest, the SGA also provides 
as special cases, exactly solvable limits of the models. 
These are known as dynamical symmetries, and can of- 
ten be found by suitable choices of the parameters of a 
model. Typically one has more forms of the Hamiltonian 
which are exactly solvable, than dynamical symmetries 
predicted from group theory M. 

The dynamical symmetry problems up to now have 
focussed essentially on the description of bound and 
scattering states. By mapping the quantum mechani- 
cal problem onto an algebraic structure, representation 
theory immediately provides the excitation energies and 
eigenstates. However, there is one missing category in 
the group theoretical description of quantum systems, 
namely, Hamiltonians with periodic potentials. In this 
case, the energy spectrum is characterized by energy 
bands and gaps. We extend the group theoretical frame- 
work by showing that dynamical symmetry and SGA 
techniques can be applied to such problems, where dis- 
persion relations, transfer matrices and energy bands can 
be understood in the context of representation theory. 



Consider Hamiltonians in one dimension of the form: 



h 2 d 2 



(i) 



For the Poschl- Teller potential, V(p) oc 1/ cosh 2 (p/p ), it 
is well known that bound states can be expressed as dis- 
crete representations of SU(2), and scattering states by 
the principal series of SU(1, 1), for which j = — 1/2 + ik, 
k > §-§. In passing from SU(2) ( where [J z , J±] = 
±J± and [J + , J_] = 2J Z ) to SU(1, 1) (obtained by chang- 
ing the sign in [J+, J—] = —2J Z ), the Casimir invari- 
ant, J 2 changes from J 2 = J 2 + J 2 + J 2 (sphere) to 
J 2 = J 2 — J 2 — J 2 (hyperboloid) , and the representation 
quantum numbers j, m fall into three well known series: 
the principal, discrete and complementary. However, if 
one has a periodic potential V(p) = V(p + a), for some 
constant a, it is not at all clear how this fits into the 
scheme provided by representation theory. In this case 
we would like to derive the transfer matrices and disper- 
sion relations. 

When one considers the eigenstates in an energy band, 
Bloch's theorem states that it must have the form 
*fe(p) = <j>k(p)exp(ikp), where (j> k (p + a) = <f> k (p) has 
the periodicity of the potential. Group theory typically 
provides single- valued wavefunctions, suitable for band 
edges. In the SU (1,1) dynamical symmetries we con- 
struct, we will see that wavefunctions can be obtained 
by (a) using multiple- valued projective representations 
(j, m) (labeled by J 2 and J z and associated with equiv- 
alence classes of vectors defined up to a phase, which is 
natural in Hilbert space), and (b) by using a basis | j, X) 
labeled by J 2 and a non-compact generator of SU(1, 1) 
such as J x . Several surprising aspects will emerge, such 
as the use of the complementary series and the need for 
non-unitary representations. 

In order to see how dynamical symmetry methods and 
representation theory might apply to Hamiltonian with 
band structure, we must first construct a dynamical sym- 
metry for a periodic potential. It is possible to do so with 
the Scarf potential 
sionless variables: 



V(p) oc 1/sin (p/po). In dimen- 



H Sc ^(x) 



A 



dx 2 



sin x 



V(x) = £$(x). (2) 



(The energy is now E = (h 2 /2Mpl)£.) This potential 
is known to have band structure for — \ < A < ||. 
By suitable choices of generators, we can realize Hs c as 
a Casimir invariant. There are different ways to do this, 



but not all are fruitful. For instance, in SU(2) we obtain 
A = to 2 -1/4, and£ = (j + 1/2) 2 . Here H Sc = J 2 + l/4is 
simply related to the Casimir invariant J 2 of SU (2). This 
realization is unsuitable since (i) the discrete representa- 
tions give only periodic solutions and (ii) the coupling A 
is never in the range of physical interest. 

If we choose the generators of ££7(1,1) in the form 
J± = e ±I *(=F sinOdg + icosOd^), and J z = —id^, we can 
realize Hs c with A = j(j + 1) and £ = to 2 : 

(3) 



Sc^jm 



TO 2 *, 



This is more natural, since a given strength of the po- 
tential corresponds to a fixed representation j, and the 
eigenvalue is related to the projection m. (We have re- 
placed 9 by x in relating Eq. || for Eq. ||) . To describe 
band spectra in SU(1, 1), we pass to the projective uni- 
tary representations which have three classes: the 
discrete series D^, with j < 0, the principal series where 
j = — 1/2 + icr, (a > 0), and the complementary series 
with — 1/2 < j < 0. While the discrete and principal se- 
ries of SU(1, 1) has found application to problems such as 
the Poschl- Teller potential discussed above, the remain- 
ing complementary series has found little or no applica- 
tion in physics. Surprisingly, this series is precisely what 
is needed for band structure, and is essential to complete 
the dynamical symmetry approach to quantum systems. 

The complementary series is shown in Fig. 1 (top) 
and corresponds to potentials with —1/4 < A < 0, 
exactly the region of interest. The unitary representa- 
tions are labeled by two quantum numbers j,m, with 
to = too ± n (n = 0,1,...), where < mo < 1 and 
mo(m(i — 1) > j(j + 1) > —i. In the figure, we plot 
a generic line j(j + 1) = const (dots) together with the 
parabola Too (too — 1) (solid). We also show the periodic 
repetition of this parabola given by (to ± n)(m ± n — 1) 
(dashes) . The values of j, to which correspond to unitary 
representations are then found to be n — j < to < j+n+l 
as well as —(j + n+l) < m < — (n — j). The relation 
of these representations to the band structure is shown 
in Fig. 1 (bottom). There we plot the eigenvalue to 2 as 
a function of j, related to the strength of the potential 
through A = j(j + 1). The range of to corresponding to 
unitary representations give the energy gaps. The non- 
unitary representations correspond to the energy bands 
(shaded region). We know that the non- unitary repre- 
sentations must correspond to the energy bands from an 
analysis of the j — > limit. Here the strength of the po- 
tential A = 0, and we have free motion, where all values 
of m must be allowed, which from Fig. l(top) corre- 
sponds to the non-unitary representations. 

The discrete (projective) representations are readily 
understood as well. The series provides the lower 
band edges, as shown by the boxes in Fig. 1 (bottom) 
for a particular choice of potential given by j (dashes). 
(Both series ± are degenerate and give the same spec- 
trum.) The upper band edges are obtained by realizing 



that the strength A = j(j + 1) of the potential is invari- 
ant under the transformation j — > — (j + 1). The dis- 
crete series D^ = 1 _ J then corresponds to the upper band 
edges (circles). The principal series has a potential with 
strength A < — i, for which the Hamiltonian is no longer 
self- adjoint, and is of no physical interest. 

The transfer matrix T can be computed directly from 
the SU(1, 1) wavefunctions (where care is taken at the 
singularity UM): 



T = 



a (3 
(3* a* 



(4) 



cos irm 



. / 2 r(^4±^)r(^p; 

in7r(i + 1) " * I m r r (-i±E) 



COSTT^COSTT^ 



f3 = ie 



2 r(^±f±i)r(^)J sin^j + i) 
/ m r(^)r(^) 



2 r(2±£2»)r( 



2+j—m \ 



1 r(^4±^)r(^p) \ cos^^costt^P 

+ - r (_irr»j r (_i±Ei) J sinTT^ + i) 

Taking cosfc7r = Re(ae mm ), we obtain the dispersion 
relation cos7rfc = cos 7rm/ sin7r(j + |), or 



£{k) 



cos 1 cos7rfc sin7r-i A 



(5) 



which agrees with that derived by Scarf ||]. The bands 
are restricted to the non-unitary parts of the complemen- 
tary series of SU(1, 1). 

Using the properties of our Scarf dynamical symme- 
try, we can consider a generalization of this potential in 
50(2,2) w 5?7(1,1) x SU(1,1). The states are labeled 
by j, m,j',c through the direct product \jm) x \j'c) of 
SU(1, 1) states. It is possible to construct generators 
of the SU (1,1) algebras, denoted J and K, for which 
j2 _ ^ g0 j — j', and which has the following 
periodic Hamiltonian as dynamical symmetry: 



gP (to + c) 2 
dx 2 



4(j + i) 2 -i 



sin 2 x 



jmc 



(Jz 



jmc 



(to - c) 2 * jr 



(6) 



(This can also be expressed as V(x) 



(a 



b cos x) I sin x.) We consider the band structure for 



•c < 



and 



< J < 



Again by symmetry, 



we use the regime — | < j < — \- To determine the 
band structure, we take the projective representations of 
50(2,2) obtained through the direct product of the re- 
sults of Fig. 1 (top). These are shown in Fig. 2, where 
we plot the complementary series for 50(2, 2) as a func- 
tion of to — c and m + c for three values of j. The shaded 
areas correspond to unitary representations. The striped 



regions are obtained from direct products of unitary rep- 
resentations, while the checkered regions are direct prod- 
ucts of non-unitary representations, which can be made 



construct the transfer matrix and the dispersion relation 
for the Lame equation with integer j: 



< 2- The representations in a = cos2k(£)K (3 = i 



unitary in the strip \m + c 

Fig. 2 correspond to (a) j — —0.45, (b) j = —0.35 and (c) 
j = —0.25. The energy gaps are indicated by the shaded 
regions. To obtain the bands for a given potential (m + c 
and j constant), one chooses j and the value of m + c in 
the range [—1/2,1/2], and draws the vertical line. The 
values of to — c which are in the unshaded regions, in- 
cluding the boundaries, correspond to the energy bands. 
One can see that as j — » — \, the bands vanish, leaving 
only discrete eigenvalues, analogous to Fig. 1 (bottom). 

The transfer matrix can also be derived (which we omit 
here for brevity), and results in the dispersion relation: 

. , . cos(2 j + lW cos(to + cW + cos(to — c)n 

cos(Trfc) = — — r4- ^- — r4 ^ — . (7) 

v ' sm(2j + l)7rsm(m + c)n v ; 

This SO (2, 2) Hamiltonian has three limits where it re- 
duces to the Scarf potential: (a) j = — |, (b) m+c = ±i, 
and (c) m + c = 2j + l. In these cases, the range of to — c 
agrees with that in Fig. 1 (with twice the period for case 
(c) ), and the dispersion relation agrees with the Scarf 
result. 

We have realized the above dynamical symmetry band 
problems using projective representations of SU(1, 1) in 
the J 2 ,J Z basis. There are other ways in which one 
can realize band structure as well, using bases defined 
through non-compact generators, such as J 2 , J x . In this 
case the states | j, A) are continuous, with A a real num- 
ber. The final case we would like to discuss here is a more 
general SGA case of this type. Consider the 1-d crystal 
formed from the superposition of Poschl- Teller potentials 
V(p\ = -(j(i + l)/Po)EnCOsh- 2 ([p-na]/p ). Through 
the linear transformation p — (iirpo/2K)x — (a + iirpo) /2, 
where K and K' are elliptic integrals, the Hamiltonian 
reduces to the Lame equation, which has a periodic po- 
tential of real period 2K : 



sm2k(£)K 



dx 2 



+ j(j + l)K 2 sn 2 (cc|Ac) 



<5 = 



(8) 



Here sn(a;|«;) =sn(x) is the doubly periodic Jacobi elliptic 
function of modulus k. This equation can be related to 
the non-dynamical symmetry problem [^|Jl^1 



H — J„ 



K 2 jy. 



(9) 



When the generators J are written in conical coordinates 
(0,4>), the above Hamiltonian decouples into two Lame 
equations, one in 6 and the other in <fi. The band edges 
are obtained from SU(2) ||, which should be contrasted 
to the previous examples where the edges were classi- 
fied according to the non-compact algebras SU(1, 1) and 
50(2, 2). For integer j, there are j + 1 energy bands. 

Non-periodic solutions of Lame's equation can be ex- 
pressed in terms of theta functions |L2] . We can explicitly 



Z)n=i sn a n dn a„ / cn a n 



(10) 



where the dispersion relation k(E) is expressed in terms 
of the Jacobi zeta function Z of modulus n: 



7T 



(11) 



The parameters a n are solutions of j coupled equations 
which depend upon the eigenvalue £ |12j : 



2 

ns ctr. 



cn a n ds a Tl 



(12) 



sn a v dn a v cn a„ — sn a„ dn a n cn a n . , . 

0=> (p^n). 

*■ — ' sn^ a r , — sn z a„ 

This transfer matrix has the following properties, (a) In 
the limit k = 1, the Hamiltonian becomes the Poschl- 
Teller potential, which is not periodic, and the denom- 
inator of (3 becomes i\J £ — j(j + 1) = ik. Here we ob- 
tain the transfer matrix with matrix elements apx = 
T(ik)T(l + ik) /r(l + ik + j)T(ik - j) and (3 PT -> 0, the 
desired result for integer j. (b) In the limit k = 0, the 
Hamiltonian (9) becomes H = J 2 while (8) is free mo- 
tion. In the basis | j, A), (9) has the desired free particle 
dispersion, £ — A 2 > 0. We also recover the free parti- 
cle transfer matrix, (c) There are additional connections 
one can make with the Mathieu equation when n — > 0, 
j — ► oo, and K 2 j(j + 1) = const., as well as connections 
to the Scarf potential. Because the equations (12) must 
generally be solved numerically, we show the type of re- 
sult one finds for j = 1. Here there is one root a, and 
the wavefunction in Bloch form is: 



H(x + a) ( mx 

^r exp \2K 



exp(— ikx), (13) 



where H and are theta functions, and a is determined 



from £ = dn 2 a + k 2 , which restricts the energy to the 



bands k 2 < £ < 1 and 1 + k 2 < £ < oo. For the lower 
band, a — K + i[3 where /3 ranges from K' to 0, while 
for the upper band a = i/3, where f3 ranges from to K 1 . 
The dispersion relation then has the form 



k(£) = -Z(p\l- K 2 ) + ^(l-J^) -£'. 



K 1 



(14) 



where £' = — £)(£ — « 2 )/(l + n 2 — £) for the lower 
band, and £' = ^J{£ - k 2 ){£ - 1 - k 2 )/{£ - 1) for the 
upper band. In contrast to the Scarf problem, we see 
here that the band edges are given by the discrete repre- 
sentations of SU(2). The bands correspond to a diago- 
nalization in the SU(1, 1) basis, such as that defined by 



J and J x . This is evident if we examine the limit k = 0, 
where we expect a continuous spectrum. 

We have shown that SGA and dynamical symmetry 
techniques can be applied to Hamiltonians with periodic 
potentials, and band structure can arise naturally from 
representation theory. This fills a long-standing gap in 
the algebraic approach to quantum systems. Hamilto- 
nians such as the Scarf's are now reduced to problems 
as simple as solving H = J%, using representations with 
multi-valued wavefunctions. The complementary series 
now finds natural application in the band structure prob- 
lem, where its non-unitary representations are impor- 
tant. Further, both compact and non-compact Lie al- 
gebras play a role in the band edges - SU(1, 1) for the 
Scarf and extended Scarf potentials, and SU(2) for the 
Lame equation. We expect that dynamical symmetries 
exist for potentials in higher dimensions, where the point 
group of the lattice will play a role in the representations 
of algebras such as SU(n,m). Within these representa- 
tions one should find the complete spectral solutions to 
the problems. 

We would like to thank F. lachello for many useful 
discussions. This work was supported by DOE grant DE- 
FG02-91ER40608. 
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FIG. 1. (Top) Complementary series of SU(1,1) labeled 
by j, m. The unitary (projective) representations sat- 
isfy mo(mo — 1) > j(j + 1), where m = mo ± n 
(n — 0, 1, ...) and < mo < 1, and coincides with the energy 
gaps. (Bottom) The corresponding Scarf energy spectrum 
E = [Ti 2 1 (2M po)]m 2 is plotted as a function of j. The band 
edges are given by the discrete series and the bands (shaded) 
by the non-unitary part of the complementary series. 
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FIG. 2. Unitary representations (shaded regions) of 
SO(2, 2) obtained from Fig. 1 (top), restricted to the physical 
region \m + c\ < 1/2, for selected values (a) j = —0.45, (b) 
j = —0.35 and (c) j — —0.25. The physical spectrum, given 
by (m — c) 2 , is obtained by choosing a value of m + c and j 
(i.e. a specific potential). The bands are the values of m — c 
which are non-unitary (in the unshaded regions). 



